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Abstract 

A one-dimensional Klein-Gordon problem, which is a physical model for a quantum 
particle submitted to a potential barrier, is studied numerically : using a variational 
formulation and a Newmark numerical method, we compute the mean position and 
standard deviation of the particle as well as their time evolution. 
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1 Introduction 

It has been well-known for a few years now that in quantum mechanics a particle can 
climb up a step even if it has not enough energy a priori and it will be reflected then 
with a delay. In classical mechanics it would just try and go back to its position. 
Everything happens here as if the particle could go through a wall (cf. Fig. 5.1 in 
|17j ) ! This phenomenon is called tunnel effect and has been a subject of interest 
for physicians and mathematicians : the delay has been measured by the physicists 
A. Haibel and G. Nimtz (cf. [IE]). It seems that it is approximately equal to the 
reciprocal of their frequency and that it is independent of the shape of the barrier. 
The Workshop on Superluminal Velocities gathered some recent contributions in 
Cologne in 1998 to describe transient phenomena concerning the tunnel effect (cf. 
|20j). In particular J. M. Deutch and F. E. Low were interested in 1993 (see |13j ) 
in superluminal effects of low-frequency Gaussian wave packets traversing a rectan- 
gular potential barrier. 

T. Hartman calculated in (TH] the tunnelling time of Gaussian wave packets for one- 
dimensional rectangular potential barriers based on the time-dependent Schrodinger 
equation. An analytic expression is given in terms of the wave number. It is easy 
to adapt his calculations to obtain an expression for the delay of reflection. 
The main obstacle for a rigorous formulation comes from the fact that T. Hartman 
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calculates pnase siiiits lor monociiromatic waves, uieariy, sucii waves cannot exniDit 
physically measurable transient features. Instead signals having a carrier frequency 
have to be studied. 

That is why F. Ali Mehmeti and V. Regnier used wave packets with a narrow fre- 
quency band together with a solution formula inspired by J. M. Deutch and F. E. 
Low in their 2003 paper ([6j and |7j|). For simplicity and to cover the relativistic case 
we chose the one- dimensional Klein-Gordon equation (as we will do in this paper as 
well) and formulated our results in terms of the energy flow to ensure the physical 
measurability of the detected transient phenomena. We gave an analytic expression 
of the delay which is in accordance with the conjecture of Haibel and Nimtz. 
More recently (cf. [8]), we got interested in causality. It is well-known that the solu- 
tion of the classical linear wave equation with compactly supported initial condition 
and vanishing initial velocity is also compactly supported in a set depending on 
time : the support of the solution at time t is causally related to that of the initially 
given condition. Reed and Simon have shown that for a Klein-Gordon equation 
with (nonlinear) right-hand side —Xu^, causality still holds. We proved in [8] the 
same property for a one- dimensional Klein-Gordon problem but with transmission 
and with a more general repulsive nonlinear right-hand side F. We also prove the 
global existence of a solution using the repulsiveness of F. In the particular case 
F{u) = —Xu^, the problem is a physical model for a quantum particle submitted to 
self-interaction and to a potential step. 

Note that this problem is a transmission problem on a simple network : two semi- 
infinite branches connected at one point. Recent results on partial differential equa- 
tions on networks, ramified spaces and more generally on multistructures, can be 
found in [U]. 

F. Ali Mehmeti solved explicitly the wave equation in a cross-shaped network in a 
seminar in 1982 ([2j). The result was then generalized by V. Regnier for any integer 
value of n in [5j. Analogous coefficients (analogous as our reflection and transmis- 
sion coefficients) appear in a similar way in the construction of prewavelets in a 
star-shaped network in D. Mercier's thesis ([23])- 

In this paper, the geometry is simpler (only two branches) but the waves are dis- 
persive. Our aim is to study the time evolution of the particle especially its mean 
position and the variance of its position. The idea has been suggested by F. Ah 
Mehmeti to V. Regnier who wants to thank him for doing so. 

In quantum mechanics neither the position of a particle nor its momentum is known 
exactly as it would be the case in classical mechanics. The Klein-Gordon equation 
which we will use in the following is satisfied by a complex-valued function such that 
the square of the modulus of this function is the probability density of the particle. 
So the mean position and its variance can give us information on the behaviour 
of the particle. Since the particle is submitted to a potential step, the potential 
energy is not constant but only piecewise constant and the problem is the linear 
transmission following one. 

We call ^1 = (— oo; 0), Q2 = (0; +00), I = [0; T] and denote by c and {k = 1; 2) 
three strictly positive constants with 02 > ai. We consider the linear transmission 
problem (P) of finding Uk '■ I x ^ IR with the time variable t G / and the space 
variable x & Qk, k = 1;2 satisfying : 
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{Ek) : ^(t, x) - c'-g^it, x) + akit, x)uk{t, x) = 0, V(t, x) E IR+* x ilk 
(To) : Mi(t,0) = M,(t,0)f Vt G G {1;2} x {1;2} 

m): ^(M-)^^(t,0-),VtG/ 
(JCi): Wfc(0,-) = /b,,A; = l;2 

(/C2): ^(0,-) = ^7b,,A; = l;2 
where / and g are defined from IR to ZR. 

Tlie equations (Ek) are Klein-Gordon equations i.e. dispersive wave equations. 
Tlie function / represents tlie initial condition and g contains the initial veloc- 
ity. This problem is analyzed in [7j. In particular the meaning of the coefficients is 
explained : the propagation velocity of the waves (phase velocity) is c (in the case 
Ofe = 0), and the Cfc's are coefficients which characterize the dispersion : the shape 
of the waves changes when time goes by. Then the velocity of wave propagation is 
not constant but strongly depends upon the frequency (dual variable in the Fourier 
transform of the signal). Recall that with classical wave equations (a^ = 0) the 
waves are only translated since the propagation velocity is constant (equal to c). 
When ai(t,x) = ai and a2(t,x) = 02, with 02 > ai > 0, the model represents a 
particle submitted to a potential step with height (02 — ai). 

(To) and (Ti) express the absence of energy loss at the central node x = (cf. [2], 
[3], [1] and Theorem 3, Section 3 of [7]). See Section 5 of the same paper for a 
physical justification and some more details about the model. 

A classical way to treat this problem which is the one we used in [7] is to use 
the spectral representation of the spatial operator A involved in the problem and 
defined as follows : 

D{A) ={ve ULiH\nk)/v satisfies (To), (T,)} 

A:ve D{A) ^ Av := ( -c^^ + akV,] in H := ULi L^^k) 



dx'^ J 



fee{i;2} 



and given u{t) = {ui{t, ■),U2{t, ■)), the above initial boundary value problem can be 
rewritten as the following abstract problem denoted by [AP) : 



(AP) 



u G C^{IR, H), u{t) G D{A), ytelR-^ 
f^{t) + Au{t) = in H 

n(0) = $,^(0) = M/ 



where $ = (/i,/2) G D{A) with : fi^ ^ IR for k G {1;2} and ^ = (c/1,^2) G 
D{y/A) with gk-Mk^ IR for k G {1; 2}. . 

The spectral theory is studied in [3], in Theorem 1.5.1. Here dj = 1, V j G {1 ■ ■ - n} 
and for us the c/s and a/s have to be interchanged. The operator A is proved 
to be self-adjoint and the existence and uniqueness of a solution of the problem 
in C2(/i?+, riLi L^^k)) n C\IR+, D{VA)) n C°(/i?+, D{A)) is stated in [3| (in the 
context of the analysis of the asymptotic time behaviour of the solution). Recall 
that D{VA) := {v G nLi^^(^fc)/^ satisfies (To)}. 
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The solution in the case gi = g2 = is recalled in the following theorem (cf. [7]). 



Theorem 1 

Assume that fk-.Q^.^ IR with k G {1; 2}, initial data of Problem (P) are such that 

{f,,h)eD{A) . 

Then the restriction to Qi of the unique solution of Problem (P) (such that {ui,U2) 
belongs to C^{IR-^ , ULi L^i^k)) n C\IR^ , D{VA)) n C^{IR+ , D (A))) is given, for 
{t, x) in IR^* X Qi, by : 



u,{t,x) = 2^/f^^_,^)CosM)^>(^(/o+°°/i(«)e-^^(-^)(— 
Recall that, for real u and j G {1; 2} 



> aj — uj 



if uj"^ < a, 



if u"^ > a, 



Note that the solution is the sum of three terms : the original term, the reflected 
one which appears due to the discontinuity in the potential and the last one, called 
"transmitted term", which contains the information on the transmission of the sig- 
nal from the second branch to the first one. 

It is the expansion in generalized eigenfunctions of cos(v^i)'^' which is the the- 
oretic solution of (AP) when vanishes. Otherwise the solution has an additional 
sine term : it is cos(A/At)$ + (v^)~^ sin(v^t)\E'. 

To study the transient behaviour of the particle we need to compute its mean po- 
sition M{t) and the variance of its position V{t) defined at the time t respectively 
by : 



M(t) 



V{t) 



Jjj^\u{t, x)\'^dx \JlR 



^x\u{t, x)\'^dx 



1 



(x - M{t)f \u{t,x)\'^dx 



Jjj^\u{t, x)\'^dx \JlR 

Using the spectral exact solution of the problem is too complicated. A numerical 
approach of the problem is thus used in this paper. 

In Section 2, the numerical scheme based on Newmark method is explained and 
stability and error estimates are given using the properties of the exact spectral 
solution. 

Section 3 is concerned with the physical interpretation of the results obtained with an 
initial condition which is sufficiently localized in space and frequency. Heisenberg's 
principle keeps the position and the momentum of a particle from being known with 
the same precision : the less uncertain is one of them, the more uncertain is the 



4 



otner one. so only compromises oasea on (oraussian lunctions can oe cnosen ana 
that is our choice. 

Our main result is that the time series of the mean position of the particle has 
a linear trend and that of the variance has an exponential one. Except from their 
sinusoidal variations (linked to the wave behaviour of the particle), the mean position 
grows linearly with time, at least while the particle has not reached the exterior node 
(reflections will then take place) and the variance (i.e. the uncertainty of the position 
of the particle) grows in an exponential way. 

2 Numerical scheme 

2.1 Variational formulation and Newmark method 

The numerical treatment of the problem, which is developed in the following, is an 
adaptation of Brezis's work about second-order evolution problems (cf. Chapter 8 
of [10], p. 197-204). 

The variational formulation of Problem (P) is, for t G / = [0; T] : 

' df{u{t),v)H + a{u{t),v) = 0, y V eV 
< w(0) = / ifeV) 
^ UtiO) =g (geV) 

where the scalar product {■', ■)h is the classical one in the Hilbert space H = 
L'^{{—L; L)), V = H^{{—L] L)) and the bilinear form a is defined, for {u, v) & V xV, 




Now the idea is to construct a sequel of spaces Sq C Si C S2 ■ ■ ■ Sm C V and to find 
a solution Um G Sm of Problem (P^) '■ 

df{Um(t),v)H + a{Um(t), v)h = 0, V f G 5'm 

(0) = gm {gm e S^) 

Since fm and gm are chosen to be approximations of / and g (discretization in the 
space variable x), Um is an approximation of u. 

Now a discretization in the time variable t is done : At = T/umax and t„ = nAt 
with n = 0; 1; 2 ■ ■ -nmax- The choices of Umax and m depend on the wanted precision 
of calculus. The numerical scheme is due to Newmark. If U"^ is Umitn) and V"^ = 
Um,t{tn) then 



(yn+i _ y)^ + (At) . a (7?7"+i + (1 - 7)[/"; v) = 

U'> = fm 

[V^ = gm 



(3 and 7 are parameters on which conditions will be put further for the scheme to 
be stable (cf. the following section.) 
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iNOW we cnoose a oasis {^j fj^Q m tne space Dm- -Lne (^ram matrix oi tnis oasis 
is Gm = (^j; the matrix of the bihnear form a is = j'- 

Denote by C", D", and the vectors containing the coordinates of t/", l^*^, 
/m and in the basis {'^Y}T=0' ^^^^ scheme is : 

' G'm(C"+i - C" - (At)/}") + (At)2A 

^ Gm(/}"+i - /}") + (At)Am (7C"+i + (1 - 7)C") = 
C = CI 

yD^ = Di 

Since it holds Um{tn,x) = Yl'j'=QC'^ipj{x), the coefficients C" contain the approxi- 
mated values of M at X = Xj, at the time t„ = nT/rimax when the ^/'/s are defined 
by ipj{xk) = Sjk for (j, A;) G {0; 1; 2; ■ ■ ■ ; m}^ with = 1 if j = A; and otherwise. 

2.2 Stability and error estimates 

The choice of adapted values for At, 7 and j3 brings stability to the Newmark 
method. 

Proposition 1 

Suppose that f G D(A^^'^) and g E D{A), /5 = 1/4 and 7 = 1/2. Then there exists a 
unique solution u E C'^^-'([0; T); D^A^^"^)) for j = 0; 1; 2; 3 and the Newmark method 
is unconditionally stable and of order 2. Denoting by Um the projection on the vector 
space Sm, the error is 

+ |(/ - nm)n(t„)| + Jo*" (|(/ - Il^)uu{s)\ + At\ui'\s)\) ds} 

The operator A and its domain D{A) are defined in the introduction. 
Proof. The first part of Theorem 8.6.2 of |10j says that, for regular enough u, that 
is, for u E C^{[0,T];V) H C'^ {[0 , T]; H) , the method is unconditionally stable (i.e. 
without any condition on At) and of second order if /5 = 1/4 and 7 = 1/2. Now 
Lemma 1.1.6 of ^ states that, if / G /}(A('=+i)/2) g^^d g E D{A''/'^), then the solution 
u belongs to C^+^-^([0; T); D{A^/^) for j = 0; 1; 2; • • • A; + 1. Since V = D{A^/'^) and 
H = D{A°), f E D{A^^'^) and g E D{A) are enough conditions. ■ 

3 Numerical results and physical interpretation 
3.1 Choice of the initial conditions 

For all the numerical computations, the velocity of the signal c is chosen to be equal 
to 1 and the potentials ai{t,x) and a2{t,x) are constant functions. On the first 
branch we choose a vanishing potential ai = 0. 

The initial condition / is a Gaussian function multiplied by the square function. This 
choice is justified by the fact that the initial condition is required to be sufficiently 
localized in space and frequency i.e. the function and its Fourier transform both 
need to be localized enough. Now the Fourier transform of a Gaussian function with 
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n+l 
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mean m ana stanaara aeviation a is a (oraussian lunction witn mean m ana stanaara 
deviation 1/a. That is why Gaussian functions are good candidates for us. Hardy's 
note about a remark of Wiener even states that they are the best candidates since 
the only pair of functions / and g (where g is the Fourier transform of /) such 
that both are 0(|x|™'exp(— 0.5a:^)) is given by finite hnear combinations of Hermite 
functions for / and g (see p!8]). 

The multiphcation by the square function then makes the function / vanish at zero. 
Thus the choice g = —cf for the initial velocity {ut{0)) makes the signal translate 
to the right until it bumps into the potential step at x = 0. Let us recall that the 
solution of the wave equation with initial conditions u{t, 0) = uq and Ut{t, 0) = vq is 

1 1 rx+ct 

- [uq{x + ct) + Mo(x - ct)] + — / vo{y)dy 

Z Ac Jx-ct 

It is the superposition of a left-going term u~{t,x) and a right-going one u'^it^x) 
defined by 

1 1 fx—ct 

W^i^, = i;Uo{x -ct) - — / voiy)dy 

2 Ic Jo 

2 \ nx+ct 

u~{t,x) = -uo{x + ci) - — / vo{y)dy 
I Zc Jo 

So with vq = —cUq and uq such that mo(0) = 0, u~ vanishes and u~^{t, x) = uo{x — ct). 
Physically an adapted initial velocity keeps the signal from "falling" and separating 
itself into two parts going to opposite directions. 

For the computations, / is the product of the square function by the Gaussian 
function with mean —3 and standard deviation 1. A constant factor is added so 
that / represents a probability density i.e. the integral over IR of / is equal to 1. 
Then it is ^ ^ 

f(x) = exp I — (x + 3)^ 

^ ^ lOv^ 2^ ^ 

and Ut{0, x) is 



9{x) = ^Q^^ x{x^ + 3x-2) exp (^-^(x + 3)^ 

The length of the branch L is chosen to be equal to 60 which avoids the reflections at 
the exterior nodes to take place too fast. They would interfere with what happens at 
the discontinuity x = 0, which is our point here. Further since the value of / on Q2 
is negligible, the probability for the particle to be on the right branch is negligible. 
So the particle is on the left branch initially (around —3.5 where its probability of 
presence reaches its maximum) and it meets a potential step. 
Note that / and g are both rapidly decreasing C°° functions. 

Proposition 2 

The functions f and g defined from IR to IR by : 

f(x) = 1= x^ exp ( — (x + 3)^^ 

gix) = ^Q^/^ x(x^ + 3x-2) exp(^--{x + 3)^^ 

satisfy f G D{A^/^) and g e D{A). 
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JrTooi. iiieorem i.o.z oi [4j gives tne equivalent statements : 

(i) u e 

(ii) u satisfies (Kq), (Ki), {K^), where 

u satisfies (Kj) u G JlLi H^+^{Vti) and 

A^^'^u satisfies (Tq), if j is even 

satisfies (Tq), if j is odd 

Tfien u G D{A) m G JlLi H^i^i) and m satisfies (To) and (Ti). Tfiose conditions 
of continuity and continuity of tlie first derivative at x = 0, are clearly satisfied by 
tlie function g wliicli is a C°° function. 

Now u G <^=^ u G nLi « satisfies (Tq) and (Ti) and Au satisfies 

(Tq). Tlie first two conditions are satisfied once more. As for the third one, it gives 
the jump of the second derivative at i.e. 

«"(o+)-«"(o-) = ^«(o) 

Since the second derivative of the function / is continuous (/ is a C°° function) and 
since /(O) = 0, this transmission condition is satisfied by /. ■ 

3.2 Interpretation of the different cases 

In quantum mechanics, the position of a particle is never known with precision. 
Only its probability density is known and given by the expression u{t, ■) at the time 
t, solution of our problem (P). We calculate the mean M{t) and the standard 
deviation a{t) of the signal u{t, ■) at the time t. Let us recall that they are defined 
as 

< 

The variations of M and a with respect to t are studied numerically here in different 
situations. 

3.2.1 Necessity of a numerical approach 

In fact an analytic approach could be envisaged in the absence of potential step 
since the situation would be less complicated. 

With a potential step (ai = and 02 > 0), and in the particular case gi = g2 = 0, 
the solution could be rewritten in terms of Fourier transforms as follows (cf. [8]) : 

Theorem 2 

Assume that fk'-^k^ IR with A; G {1; 2}, initial data of Problem (P), described in 
the introduction, are such that (/i, /2) G D{A), fi is compactly supported in (—00; 0) 
and /2 = 0. 

Then the restriction to Vti of the unique solution of Problem (P) (such that {ui,U2) 
belongs to C\IR^ ,Ul=i L\nk)) n C\IR-^ , D{VA)) n C^{IR+ , D{A))) is given, for 
{t, x) in IR^* X Qi, by : 
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Ui{t,x) = —J^^^^ cos("\/ c^uj^t)Tfi {uo) 



IT 



where the complex square root has been defined such that ^ — 



i\Ja2 — (P'LO^ if — 



— vAu^ 

Likewise, for (t,x) in IR'^* x Q2, U2{t,x 



^ <UJ< 

02 if uj<-.i^ 



271 J- 



cos 



{V(Kj^t)t 



2UJC 



Tfi {uS) duj 



So analysing the time evolution of the mean and of the variance is not simple. If 



the potential is constant {a\ = a2 

1 



the solution is : 



ui{t,x) = -J^.^^ 



TT 



with the initial condition /. Rewriting the cosine as half the sum of two complex 
conjugated exponentials leads to : 



ui{t,x) 



where uj{^) = ^/m^~+~c^. The expression {^x — uj{C,)t) is the phase of the first 
integral and it is stationary when x/t = 00' (C,)- Moreover 



Thus u is an increasing function. If the initial condition is chosen to lie in the 
frequency band [^i;^2] i-e. if the support of J-'f is a subset of [^i;^2], then at a 
fixed time t, the solution u{.,t) lies essentially in [a;'(,^i)t; a;'(^2)^]- Since the center 
and the amplitude of this interval both grow linearly with t, the mean value and 
the standard deviation of the particle submitted to a constant potential are due to 
grow linearly with time. This approach is useful to study the time decay rate of the 
solution of the Klein-Gordon transmission problem. Marshall, Strauss and Wainger 
were interested in that problem (cf. [2^), F. Ali Mehmeti as well (P] and ^). 



Note that the phase of the second integral can be obtained from the first one, 
changing t into — t. It can be seen as a negative-energy particle solution term prop- 
agating backwards in time which describes a positive-energy antiparticle solution 
term propagating forwards in time. This reinterpretation is due to Feynman and 
explained in [Ij. In fact negative-energy solutions to the Klein-Gordon equation 
may appear due to its construction (cf. Section 6.2.4 in [?]) and Feynman has given 
them a physical interpretation. 

In our situation, the potential is not constant (02 — ai > 0) so the time evolution of 
the mean position and of its standard deviation is not that clear. This justifies the 
use of a numerical scheme to study the problem. 
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d.z.z ivietnoaoiogy 

The numerical scheme, due to Newmark and recalled previously in Section 2.1, has 
been implemented in a C~^^ program. The conjugate gradient method has been used 
to solve the system. In each case the graph of the solution u(t,-) is produced for 
different values of t, as well as that of the mean position M and standard deviation 
a (as functions ot the time t). 

Two aspects are interesting : 

• first the behaviour in the neighbourhood of x = 0. 

How does the particle interact with the barrier it meets ? Different situations 
have to be looked at thoroughly : since the initial condition is fixed here, 
the energy of the particle is fixed and depending on the height of the barrier, 
whether it will have enough energy to climb it up in a classical way (as a ball 
would do meeting an inclined plane in classical mechanics) or it has not, but 
will go through the barrier due to the specific laws of quantum mechanics. 
The limit case is of course particularly interesting. 

Let us recall how the energy of the particle is defined. The initial condition 
/ is the probability density of the particle and the Fourier transform of / 
gives the frequency distribution of the signal. The energy of the particle is 
proportional to its frequency when it is sufficiently localized in frequency. The 
Fourier transform of / is / defined by : 



where Const ^ 1.111 ■ 10~^. Its shape is that of a mexican hat centered at —3. 
So the particle has enough energy to climb up a barrier with height inferior to 
9 (the square of 3) and the tunnel effect takes place for higher barriers. 

• second, the behaviour for large values of t. 

In fact not that large since we have to stop the experiment just before the 
wave hits the boundaries to avoid additional reflections. Numerically the ini- 
tial condition has been chosen far enough from the boundary so that this does 
not happen too quickly. 

Since the mean value and the standard deviation may both have a wave be- 
haviour (the solution has a wave behaviour), we study the trend of their evo- 
lution using Dagnelie's kernel smoothing (cf. Section 15.4.2 of [H]). The idea 
is to keep only the trend of the phenomena removing their local variations as 
it is done with time series. The effects of seasonality are removed to be able 
to see the evolution from one year to the following one. 

Our conjecture is that, for all the shapes of barriers (infinite, finite, piecewise 
constant or not) we envisage, both trends are linear. That is checked here using 
a computation of the correlation coefficient between the mean (respectively the 
standard deviation) and the time. The equations of the least squares straight 
lines are also calculated. The notations are : mean = At + B and r is the 
correlation coefficient of this regression, and standard deviation = Ait + Bi 
where ri is the correlation coefficient (cf. Table 1 of Section 3.2.3). 
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^.z.a iiesuits ana pnysicai interpretation 

The observation of the different graphs of the solution, its mean value and its stan- 
dard deviation leads to general remarks about the common behaviour of the solution 
(independently of the shape of the barrier) : 

The standard deviation is always constant till the particle hits the barrier and the 
mean position grows in a perfect linear way. This is due to the absence of dispersion 
on the left branch : there is a simple translation of the signal without any reshaping. 
We observe that both the mean position and the standard deviation have a chaotic 
behaviour when hitting a high barrier, which is of course not surprising. They are 
both smoother after the shock. 

In particular the standard deviation decreases a bit before increasing again and 
keeps increasing then (see Figure 1). Keeping the analogy with a ball being thrown 
against a wall, we could say that this small decreasing in the standard deviation can 
be seen as a reshaping of the particle which is somehow crushed, a bit flattened like 
a rugby ball when it hits the barrier, exactly like a round ball, which would become 
imperceptibly oval due to the shock. 

The analogy is not perfect however : the decreasing of the standard deviation is 
not larger for particles having more initial energy (as could be expected). It is ex- 
actly the same one for a fixed height of the barrier. In classical mechanics a ball 
thrown with higher velocity (and so more kinetic energy) will get more flattened 
when hitting the wall. In fact the influent parameter here is the height of the bar- 
rier : the standard deviation is larger for higher barriers. What happens is linked 
to the group velocity of the signal. On the left branch it is c, and on the second 
branch it is given by (1) with m = a2- Thus an increasing height of the barrier 02 
implies a smaller velocity of the particle on the left branch. So the signal on the 
left branch keeps moving to the right (reducing its standard deviation, but indepen- 
dently of the height of the barrier) while the right part of the signal moves slowly 
and so increases its standard deviation slowly (and the higher the barrier the more 
slowly it increases). The combination of both phenomena makes the decreasing in 
the standard deviation bigger for higher barriers. Everything happens as if a higher 
barrier corresponded to a "tougher" wall. 

Another way to explain that the decreasing of the standard deviation is not larger 
for particles having more initial energy is to say that if the height of the barrier 
is fixed while the initial energy of the particle increases, then the transmitted part 
increases. A part of the energy is used by the particle to cross the barrier partially 
and so it is not flattened as much as if it was completely stopped by the barrier. 
Note that in all cases, the standard deviation is affected by the barrier sooner than 
the mean position. The amplitude of its variations is also larger. In fact the mean 
keeps on increasing (but a bit more slowly) while the standard deviation decreases. 
The shock seems to have a small impact on the mean position : a peak during less 
than 0.05 units of time while the standard deviation is perturbed for about 0.3 units 
of time. 

To finish with the general results, observed for any value of a2. Table 1 gives the 
correlation coefficient, denoted by r (respectively ri), between the mean M (re- 
spectively the standard deviation a) and the time t (cf. Section 3.2.2). As it was 
conjectured in the introduction, both have a linear behaviour for large values of t. 
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it was not tnat clear a prion lor small Darners. 
Let us now consider the different cases : 

• Infinite barrier with small height 02 < 9 : 

In that case, the energy of the particle is enough for it to climb up the barrier 
in the classical way and the smaller the barrier, the easier it is for it to climb 
up. So the mean position which starts at about —3.5 (the particle is on the 
left branch at the beginning), keeps increasing and takes positive values after 
some value of t (denoted by to) which depends on the height of the barrier. 
The higher the barrier is, the larger to becomes. In fact the higher the barrier, 
the larger the reflected part is (cf. the values of the reflection coefficient repre- 
sented in Figure 2 of [7j) and if the reflected part is high, it means that a part 
of the signal is on the left branch, which makes the mean position smaller. 
That is why to is higher. 

The height of the barrier influences the behaviour of the particle when hitting 
the barrier : the peak of the mean position is smoother for a smaller barrier 
and the increasing of the propagation velocity is not so big for a small barrier 
since the peak is larger. Note that this apparent acceleration of the particle is 
in contradiction with the decreasing of the propagation of the velocity on the 
right branch compared to the left one (cf. (1) once more). In fact this accel- 
eration of the mean position takes place exactly when the standard deviation 
increase in a brutal way so the mean is less representative of the reality. The 
uncertainty on the position of the particle increases brutally. 

Also note the wave behaviour of the standard deviation after the shock while 
the mean position decreases perfectly linearly (except for very small barriers 
where small oscillations take place). The standard deviation increases more 
quickly for low barriers than for high ones, reaching at the same time (15 units 
of time) about 6.2 if 02 = 2, about 2.75 if a2 = 2.5, about 1.56 if a2 = 4, about 
1.05 if 02 = 5, about 0.71 if 02 = 6. Note that the oscillations are more and 
more numerous (from 5 if 02 = 2 to 7 if 02 = 6) but less and less high. The 
signal is quite flat for 02 = 1 and a2 = 8. 

If the height of the barrier tends to zero, the situation tends to that of the 
classical wave equation (no dispersion) so the signal tends to translate without 
changing its shape. 

• Infinite barrier with large height 02 > 9 : 

In that case, the energy of the particle is not sufficient for it to climb up the 
barrier in the classical way. The tunnel effect takes place to allow the particle 
to go through the barrier. A previous energetic study has shown that in fact, 
the particle lingers on in the neighborhood of x = for some time, called 
delay of reflection, before it is reflected (cf. [7]). We observe that, except for 
a delay, the energy just after the shock is equal to what it was just before. 
That delay is visualized by a spatial shift (called Goos-Hanchen shift) in the 
two-dimensional case of the double-prism experiment of Haibel and Nimtz (cf. 
[16J and [2S]). 

After the chaotic behaviour due to the shock with the barrier, the mean po- 
sition decreases regularly and linearly and the symmetry of the curve means 
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tnat tne renection is total exactly as we nave just expiamea [lae slope oi tne 
curve after the shock is opposite to that of the curve before the shock). 
Furthermore, the more energy the particle has (compared to the barrier it 
meets), so the lower the barrier, the smaller the delay is (cf. Figures 3 and 
5 of [7J). When 02 increases, the propagation velocity on the left branch de- 
creases so we have conjectured that the particle lingers on for a longer time 
around zero before being reflected and that it is the reason why the delay is 
bigger. 

Observing the shape of the mean position curve for 02 varying between 9 and 
100, we detect a peak which is larger and larger as 02 increases. About 0.1 
units of time if 02 = 9, about 0.25 if a2 = 15 and 0.3 for 02 = 100. For a2 = 9 
the mean position increases with velocity c (i.e. 1 in our numerical computa- 
tions) till the wave front reaches the potential step. Then it keeps increasing 
but with higher velocity and decreases with the same velocity after a very 
small time before decreasing with velocity c again. Now if 02 increases, the 
mean begins with decreasing (from —0.7 to —0.72 for 02 = 9.2, from —0.79 to 
—1.52 for a2 = 12 and from —1 to —1.7 for 02 = 100), before having the same 
behaviour as for 02 = 9, except for the maximal value which decreases with 
increasing 02 (about —0.55 if 02 = 15 and —0.8 if 02 = 5000). The transmitted 
part is smaller so the mean is smaller as well. 

This observation confirms our conjecture : the particle tries to go through the 
wall and the interactions with the barrier make it linger on for a while, hence 
the delay in the reflection. 

As for the standard deviation, it has the same behaviour before the shock 
as what it was for small barriers but afterwards, it becomes constant again 
with the same value as for t = 0. This is due to the total reflection and the 
absence of dispersion (the potential vanishes) on the left branch. 

• Other shapes of barrier : 



Table 1 (Potential step) 





a2 = 2 


0-2 = 6 


0-2 = 9 


02 = 15 


a2 = 150 


A 


-0.9114 


-0.9928 


-0.9940 


-0.9941 


-1 


B 


-2.4414 


-3.5840 


-3.7679 


-3.9338 


-4.2629 


r 


-0.9999 


-0.9999 


-0.9999 


-0.9999 


-1 


Ai 


0.3266 


0.0212 


0.0013 


2.4873 ■ 10-6 


-2.23 ■ 10-6 


Bi 


1.5078 


0.6682 


0.6533 


0.6525 


0.653 


ri 


0.9984 


0.9945 


0.9905 


0.9820 


-0.998 
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Standard deviation versus time for a potential barrier with height 10. 
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Figure 1: commentaires?? 
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Standard deviation versus time for a potential barrier with height 150. 
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Mean position versus time for a potential barrier with height 10. 
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Mean position versus time for a potential barrier with height 150. 
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Standard deviation versus time for a potential barrier with height 6. 
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